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Fusion mechanism in fullerene-fullerene collisions 
—The deciding role of giant oblate-prolate motion 
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We provide answers to long-lasting questions in the puzzling behavior of fullerene-fullerene fusion: 
Why are the fusion barriers so exceptionally high and the fusion cross sections so extremely small? 
An ab initio nonadiabatic quantum molecular dynamics (NA-QMD) analysis of C 60 +C 60 collisions 
reveals that the dominant excitation of an exceptionally “giant” oblate-prolate H g ( 1) mode plays the 
key role in answering both questions. From these microscopic calculations, a macroscopic collision 
model is derived, which reproduces the NA-QMD results. Moreover, it predicts analytically fusion 
barriers for different fullerene-fullerene combinations in excellent agreement with experiments. 

PACS numbers: 36.40.-c 


INTRODUCTION 


After the pioneering Cgo+Cgo collision experiment of 
E.E.B. Campbell et.al. in 1993 Q, cluster-cluster colli¬ 
sions became a versatile new field of research (for reviews 
see with lasting interest @-0|. In particular, fu¬ 

sion between fullerenes has been studied in great detail, 
both experimentally and theoretically [l, 21- 3(3] • Fusion 
is a universal phenomenon in collisions between complex 
particles covering many orders in size and energy from 
heavy ion collisions in nuclear physics 38[ to macroscopic 


liquid droplets El El or even colliding galaxies [ilj] . It 
is a great challenge of ongoing interest to reveal universal 
similarities and basic differences of these mechanisms. 

Usually, the gross features of fusion can be understood 
with macroscopic arguments 0, EM El liM] leading to 
the general expression for the fusion cross section a as 
function of the center mass energy U c . m . = E of 


one order of magnitude larger (around 80 eV) [29] and 
the cross sections are even two magnitudes smaller (a few 

A 2 ) 0. 

Up to now, there is no definite explanation for these 
findings, albeit some possible phenomenological reasons 
have been discussed 0-0, Ell, EM] • In addition, previous 
(at that time still approximate) microscopic Quantum 
Molecular Dynamics (QMD) calculations predicted the 
large fusion barriers [28j. From these studies it is also 
well known that only very few mutual initial orientations 
of the colliding cages lead to fusion (without identifying 
them). Anyway, fullerenes typically do not fuse if they 
touch, even at high impact energies, and the question 
remains, why? In this work, we provide a clear answer 
to this longstanding question. 


MICROSCOPIC RESULTS 


a(E) = 7 tR 2 


1 - 


Vb 

E 


(1) 


with R 12 = Ri + R 2 the sum of the radii of the collid¬ 
ing partners and Vb the fusion barrier (for a derivation 
see e.g. 0, |38. M)- This formula (known as “critical dis¬ 
tance model” in nuclear physics [38] or “absorbing sphere 
model” in chemistry 0 ) describes quantitatively the ex¬ 
perimental fusion cross section for atomic nuclei (with 
Vb > 0) and liquid droplets (with Vn = 0) 0, ES 0 
in the low energy range with E >Vb [44J. It is expected 
to hold also for collisions between metallic clusters [id] 
(even with Vb < 0); see last paper in the series 0-0. The 
physics behind formula m is indeed simple: fusion takes 
place, if the colliding partners touch, owing to the larger 
binding energy of the fused compound |42|. For collid¬ 
ing fullerenes one naturally expects a fusion barrier Vb 
of at least the sp 2 bond breaking energy, or more perti¬ 
nently the Stone-Wales transformation energy [33j, |34j of 
a few eV [IB] and, according to ©, a fusion cross section 
of the order of the geometrical one, 7r R 2 2 ~ 150 A 2 0. 
Experimentally, however, the fusion barriers are about 


Motivated by our recent findings of the dominating 
role of the A g {\) breathing mode in Ceo-laser interac¬ 
tion [46l . |47j, we have reanalyzed fullerene-fullerene col¬ 
lisions with the help of the ab initio nonadiabatic quan¬ 
tum molecular (NA-QMD) method iM-5l|. For systems 
as large as we are investigating here, NA-QMD is nu¬ 
merically more efficient than its ab initio QMD approx¬ 
imation H 511. In extension to previous studies 0E3- 
[29l . EM] we include a normal mode analysis 521 of the 
vibrational kinetic energy. This method decomposes the 
total kinetic vibrational excitation energy E v ,b into the 
individual contributions of all 174 eigenmodes of Ceo as 
function of time t according to 


60 174 / 60 


1/=1 


with the atomic carbon mass m c , the atomic velocities 
rj (in the molecular center of mass system without ro¬ 
tational components) and the normal mode eigenvectors 

1 * 1 v ■ 
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FIG. 1. (color online) Normal mode analysis of central 
Cgo+Ceo collisions at an impact energy of E c . m . = 104 eV. 
Shown are the vibrational kinetic energies E v ib of all 174 in¬ 
ternal normal modes as function of time t for (a) typical scat¬ 
tering and (b) typical fusion events. Distances of closest ap¬ 
proach R Tet are indicated by vertical, dashed lines. The dom¬ 
inating role of the H g ( 1) mode (black thick lines) is apparent. 


In fig. [□ such an analysis, is shown for two cen¬ 
tral Cgo+Ceo collisions with the same impact energy of 
E c m = 104 eV but for different initial orientations of 
the clusters, leading in one case to scattering (fig. Oja)) 
and, in the other, to fusion (fig. [3b)). In both cases, the 
extraordinary dominance of the H g ( 1) mode is obvious. 


During approach, this mode is very strongly excited to 
a much higher degree than any other vibrational mode. 
Its excitation energy of several eV is huge as compared to 
a single quantum of this mode of hu! mo de ~ 34 meV [53] 
and, thus, its amplitude is extremely large, as compared 
to a typical displacement of the elementary excitation 
(“giant” H g ( 1) mode). Consequently at the distance of 
closest approach R let , a highly deformed oblate-oblate 
configuration of the double cluster system is formed (see 
above illustrations in fig.[l|). This clearly distinguishable 
state accommodates practically the whole impact energy 
into deformation (potential) energy, which is quantita¬ 
tively shown in fig- Eta) , also for other impact energies. 
Up to this stage of the collision, there is no appreciable 
difference to the other collision systems (nuclei, droplets), 
where at this “critical distance” the system loses immedi¬ 
ately its memory and the energy is dissipated into inter¬ 
nal degrees of freedom (DOF) leading to a hot compound. 

The fundamental difference to the other systems con¬ 
sists in the specific properties of the oblate-prolate mode 
in fullerenes and its special role in collisions. First of all, 
among all vibrational modes, the H g ( 1) mode in Cgo has 
the largest oscillation period of T = 122 fs [53]. This 
is comparable with a typical collision time and there¬ 
fore, once excited, this mode will not lose immediately 
its memory, as clearly seen in figs. [TJ a) andOJb). Sec¬ 
ond, the oblate-prolate mode is the only eigenmode which 
can couple directly to the relative motion via its elon¬ 
gated prolate-prolate configuration, provided it survives 
the dissipative coupling to all the other internal modes. 
This is exactly what happens in Cgo+Cgo collisions and 
results in the majority of cases in scattering, like a fis- 




FIG. 2. (color online) Typical NA-QMD results for central 
C 60 +C 60 collisions at impact energies of E c . m . = 40 eV (red), 
60 eV (blue), and 104 eV (black), indicated as dotted hori¬ 
zontal lines, and obtained by ensemble averaging over 20 dif¬ 
ferent initial orientations, (a): The potential (deformation) 
energy of the relative motion U(R) in the entrance channel 
as function of the distance R. The locations and lengths of 
the vertical arrows denote the distances of closest approach 
R Tet and the total deformation energy E Tet = U(R Tet ) stored 
at R Tet , respectively. The dashed line is the harmonic fit 
U(R) = ^ (R — R°y with k = 0.12 a.u. and R° = 18.5 a.u. 
(b): The kinetic energy of the relative motion Ekin as func¬ 
tion of time t for scattering events. The mean values of the 
final kinetic energies A E ~ 17 eV and of the corresponding 
energies dissipated into internal DOF Ed(oo) = E c . m . — A E 
are indicated by double arrows. 


sion process with the prolate-prolate configuration at the 
scission point [42J (note the pronounced excitation of the 
mode in fig. Oja) at t ~ 165 fs). Only strong coupling to 
the internal DOF can prevent this mechanism, allowing 
for fusion (note the strongly damped prolate oscillation 
of the mode in fig.[D]b) at t ~ 165 fs). Thus, the compet¬ 
itive coupling of the oblate-prolate mode to the relative 
motion and to all the other (bath-like) vibrational DOF 
determines the reaction channel ! 


The strength of the bath coupling is solely determined 
by the amount of energy stored in the mode at R let . This 
coupling dominates, if the energy exceeds a certain limit 
which generally can happen only at appropriate large im¬ 
pact energies. This explains (at least preliminarily and 
qualitatively) the high fusion barriers. At a fixed impact 
energy just above the barrier E cln , > Vb (as in fig. [TJ), 
only very rare and specific initial orientations of the clus¬ 
ters can lead to high H g ( 1) excitation energies , namely, 
those with the principal axes of the H g ( 1) mode aligned 
to the collision axis, ensuring maximal energy transfer, 
which is the case in the example shown in fig. |3b) (note 
that -E v ib of the mode during approach in the case of 
fusion (fig. [3b)) is twice as large as compared to scat¬ 
tering (fig. Oja))). This, finally, explains the low fusion 
cross sections and completes the present, new picture of 
the fusion mechanism. 


It modifies also the hitherto existing interpretation of 
scattering, as a “bouncing off” mechanism [j, [j, 22, 28, 
[3lj ], like in collisions between two soccer balls. Instead, 
“fission” via the prolate-prolate configuration strongly 
suggests, that the final kinetic energy of the fragments 
is largely independent on impact energy. This is nicely 
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FIG. 3. (color online) Phenomenological collision model with 
two DOF, the distance between the centers R(t) and the di¬ 
ameter D(t) of the fullerenes: Two collinear colliding springs 
(with spring constants fc mo d e and masses m) interact during 
approach via a third, massless spring (spring constant k and 
initial length _R°) located at the collision center. 


confirmed in the calculations and shown in fig. ©b). 

The collision scenario presented in figs. [T] and [2] is 
characteristic for fullerene-fullerene reactions and qual¬ 
itatively observed in our NA-QMD calculations also for 
the other combinations and finite, small impact param¬ 
eters. Despite its microscopic complexity, the mecha¬ 
nism is nevertheless simple and can be understood and 
described by ordinary macroscopic concepts, as will be 
shown in the following. 


COLLISION MODEL 


The basic idea is to reduce drastically the 360- 
dimensional scattering problem to a one-dimensional one 
with only two, but relevant collective DOF, treated ex¬ 
plicitly in a time-dependent fashion: the distance be¬ 
tween the centers R(t) (relative motion) and the diame¬ 
ters of the fullerenes D(t) (aligned along their principal 
axes of the oblate-prolate mode). Both are coupled via 
the contact distance R c = R — D (see fig. ED- The cou¬ 
pling to the other internal DOF will be treated implicitly 
in the exit channel only. The macroscopic model is de¬ 
signed as follows: 

(i) In the entrance channel, the system consists of two 
collinear colliding springs with initial lengths D(t = 0) = 
D°, spring constants fc mo de and masses m = 4r at the 
ends (describing the fullerenes with mass M and their 
H g ( 1) modes). Tightly located in between there is a 
third, massless spring with initial length R° c and constant 
fc, describing the repulsive potential U during approach 
(remember fig. ©a) and see fig. [3] and fig. Oja)). 

(ii) In the exit channel, the massless repulsive spring is 
replaced by a “dissipative” potential Ud, which describes 
the coupling to all other internal vibrational DOF, and 
hence, controls the reaction channel (see fig. ©b)). 

The coupled Newton equations in the entrance channel 


FIG. 4. (color online) Schematic plot of the potential energy 
in the entrance channel U(R C ) (left) and, for the case of scat¬ 
tering, in the exit channel U d {Rc) (right) as function of the 
contact distance R c . Dotted blue lines with directional arrows 
indicate a typical trajectory (idealized as a straight line in the 
exit channel) with impact energy E c .m. at the initial distance 
Rc and potential energy E let at the distance of closest ap¬ 
proach R™ . The asymptotic value of the potential energy in 
the exit channel Ud(R c —I oo) = Ed(oo) is also indicated. 


read 


pR = 

Pmode-D — 


dU 

dR c 


R C =R—D 


dU 


dVmode . f 
dD + 2 dR, 


( 2 ) 


R c —R—D 


with the reduced masses p = ^ and the mass con¬ 
stant p ni ode — X - The harmonic potentials are given by 
U{R C ) = | (R c — Rc) 2 with k = 0.12 a.u. (see fig. ©a)) 
and Vmode(-D) = fcm ° de (D — D 0 ) 2 with the spring con¬ 
stant fc mo de = Pmode w^ ode = 0.51 a.u., obtained from 
the experimental frequency of the H g ( 1) mode in Ceo 
(wmode = 273 cm -1 [53|). The equations of motion 
(EOM) © can be solved analytically (see appendix). 
They describe the collision up to the distance of clos¬ 
est approach, i.e., the classical returning point R T c et . At 
this point the potential U(R C ) for R c > R T c et is replaced 
by the “dissipative” one U d {R c ) which, in dependence 
on Ec.m., controls the outcome. It is therefore repulsive 
(leading in any case to scattering) or attractive (leading 
usually to fusion, but not necessarily always). Thus, it 
has the general form 


U d (Rc) = (E ret - E d { 00 )) f{R c - RT) + E d ( 00 ) (3) 

where the form factor / must fulfill the conditions, /(0) = 
1 (ensuring the continuity of the potential at R r c et ) and 
/(oo) = 0 (making sure that the maximal amount of 
dissipated energy cannot exceed E d ( 00 ), in the case of 
scattering). In fact, the concrete radial dependence of 
the potential © is not relevant, and thus, we choose a 
simple exponential form of f(x) = exp (—^) with A = 

| R ' c ^ Rc (1 — ^ ) |, which guarantees also continuity 

of the force at R r c et in the case of scattering. 

With this, all model parameters (fc, fc m ode, A E) are 
fixed and the EOM © can be easily solved numerically. 
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FIG. 5. Comparison between NA-QMD (dashed lines) and 
model calculations (solid lines) for central C 60 +C 60 collisions 
at three impact energies E c . m.: distance between the centers 
R in a.u. (top), kinetic energy of the relative motion Ekin 
(middle), and vibrational kinetic energy of the H g ( 1) mode 
E v ib (bottom) in eV as function of time t. 


The results are shown in fig. [5] and compared with NA- 
QMD calculations (for movies see www.dymol.org). 

The model reproduces nearly precisely the microscopic 
calculations for both scattering (E c _ m , = 40, 60 eV) and 
fusion (E c . m . = 150 eV). The ongoing oscillations of some 
quantities in the exit channel are the natural consequence 
of the absence of a damping mechanism for the H g ( 1) 
mode in the spring model. The most impressive result, 
however, concerns the fusion barrier predicted by the 
model of Vb = 85 eV, which is in excellent agreement 
with former (extremely expensive) QMD calculations [28] 
of V B = 80 eV. 


ANALYTICAL SOLUTION FOR THE FUSION 
BARRIER 


The simplicity of the EOM @ and the transparency 
of the ansatz 0 encouraged us to further elaborate the 
model and to demonstrate its predictive power. As men¬ 
tioned already, the EOM f2]) can be solved analytically, 
and, from this solution one can derive an approximate ex¬ 
pression for the energy stored during approach E let (see 
FigHta)), resulting in a linear dependence on the impact 
energy of E let = a(n) E c . m . with k the ratio of the spring 
constants 


k = 


^mode 


k 


M, ,2 

4 ^mode 


k 


( 4 ) 


and the universal function a(n) of 


a(n) 


1 

4 (k 2 + 1) 


K — 1 — \/ k 2 + 1 

\/ K T 1 + V K 2 +7 


X 


7 r k 
sm(-- 


+ 1 + V K 2 + 1 


) 


K — 1 + V K 2 + 1 

\J K + 1 — V K 2 +7 

( 5 ) 


(see appendix). The difference (E ret —E d ( 00)) determines 
the (positive or negative) slope of the potential U d ([3]) at 
the returning point R r c et . Taking an idealized trajectory 
as shown in fig. [I] (i.e., neglecting the acceleration of the 
mode near the barrier), the fusion condition simply reads 

E iet = E d ( 00) at E c , m. = V B - With E d ( 00) = E c , m -AE 
(see fig. 0(b)), the fusion barrier becomes 


V B 


A E 

1 — cn(re) 


( 6 ) 


For C 60 +C 60 collisions, this approximate expression gives 
Vb = 82 eV, which is very close to the exact model 
value (Vb = 85 eV) obtained in the upper dynami¬ 
cal calculations. To obtain a first insight about the 
qualitative trends for the other combinations (C 60 +C 70 , 
C 70 +C 70 ), we ignore subtleties and use the same k and 
A E values as obtained already for Ceo+Ceo by NA- 
QMD fine tuning. The internal spring constant for C 70 , 
however, is carefully chosen and fixed again by exper¬ 
iment. From the (partly) non-degenerate five H g ( 1) 
modes in C 70 ( E ' 2 , E", A j) an experimental mean value of 
w mo de(C 7 o) = 261 cm -1 has been reported [54( giving the 
spring constant k mo( i e (C 7 o) = 0.55 a.u. For the (slightly) 
asymmetric C 60 +C 70 collision a reasonable mean value 

of k m ode(Ceo) and /''mode (C7O) of /I'rriode (Ufj() /U 70 ) — 

0.53 a.u. is used. 

With these parameters, the predicted fusion barriers 
from J 6 ]) are compared with high precision QMD val¬ 
ues [28(, in table Q] (first two rows). The analytical model 
reproduces the right trend and delivers absolute values 
within 10% accuracy. Obviously, the H g ( 1) frequencies 
w mo de and fullerene masses M, eq. J4]) determine the fu¬ 
sion barriers in fullerene-fullcrene collisions. 

This is strongly supported by comparing the predic¬ 
tions with experimental data. In this case, the finite 
temperature (T « 2000 K) of the colliding fullerenes has 
to be taken into account 29]. This has been done in the 
former QMD calculations [28| and led to the (well known) 
perfect agreement with the experimental data (cf. third 
and fifth rows in table Q]). In the present model, a finite 
temperature can be naturally taken into account by re¬ 
ducing the mode frequencies w mo de- Using an arbitrary 
common scaling factor for fc mo( je of 0.85, the predicted 
fusion barriers by eq. ([ 6 ]) are in beautiful agreement with 
the experiment (fourth and fifth row in table [T]). 


CONCLUSION 

To summarize, ab initio NA-QMD studies have finally 
cleared up the fusion mechanism in fullerene-fullerene 
collisions. The excitation of a “giant” H g ( 1) mode ex¬ 
plains both large fusion barriers and small fusion cross 
sections. This microscopic picture is non-ambiguously 
confirmed by a macroscopic spring model which depicts 
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TABLE I. Fusion barriers Vb in eV for various fullerene- 
fullerene combinations as predicted by our former QMD cal¬ 
culations [HI (first and third rows) and the present analytical 
model (second and fourth rows). The experimental values [2§] 
are presented in the last row. 



T = 0 K 

T 

= 2000 K 


QMD model 

QMD model 

exp. 

C 60 TC 60 

80 

82 

60 

65 

60±1 

C 60 TC 70 

94 

87 

70 

69 

70T6.5 

C 70 +C 70 

104 

93 

75 

73 

76±4 


clearly the physics, reproduces the NA-QMD results and 
the experimental fusion barriers quantitatively. 

We note finally, that a “giant” vibrational excitation 
of the A g ( 1) breathing mode in Ceo has been found re¬ 
cently in a time-resolved laser experiment fitl ]. The gen¬ 
eral investigation of large amplitude motion in fullerenes, 
including laser-induced fission |[47{], could become an in¬ 
teresting new field of research. 


APPENDIX: ANALYTICAL SOLUTION IN THE 
ENTRANCE CHANNEL 


Inserting the harmonic potentials U(R C ) and V(D), 
the EOM ([2| can be written as 

/rl? = — k ( R c — Rc)\r c ^r-d 

T mode-D = —kmode (D — D°) + — ( R c — R (!) 

(A.l) 

which can be solved by making an exponential ansatz 
e lQt . Doing so, the fundamental eigenfrequencies fl\/ 2 of 
the system read 

fl\/ 2 = LO \]K + 1 ± \/K 2 +7 (A.2) 

v. v ✓ 

/l/2 (M 


with the frequency a; = y ^7 

k as defined in eq. 

With the initial conditions 

i?(0) = R° = D° + R° c , 
D( 0) = D°, 


and the force constant ratio 


R( 0) = V c .m. 
D( 0) = 0 


the solution of (IA.1D is given by 

R(t ) = Vc ' m ' y ai sin (flfi) + R°, (A.3) 

D(t) = Vc m ' bi sin (flfi) + D° (A.4) 

w *= 1,2 


with the amplitudes 04/2 
&i/2 = Tj 7^ y=y 


2 /l/2(^) 


(it 




+1 


= ) and 


With the analytical solution (IA.3I) . (IA.4I) the potential 
energy at the returning point E let can be calculated. E let 
is given by 


E^^U(RT) = ^(RT-R°c) 2 (A.5) 

with R T c et = R c (t iet ) = R(t iet ) - D(t iet ). 

The returning time f ret defined by R(t let ) = 0 is approxi¬ 
mated by f ret ~ The approximation is justified since 
the second term of the sum in eq. (IA.3I) dominates for the 
parameter range of k used here (^ <C 1). 

Inserting the analytical solution (IA.3I) . (IA.4I) in the defi¬ 
nition (I A. 51) and using the relation (IA.2I) we find 

E ret = a(K) E c . m . 

with the impact energy E c m = ^ and the coeffi¬ 

cient a(n) as defined before in eq. ®. 
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